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ABSTRACT 

We present high spatial resolution maps, obtained with the Plateau de Bure Interfer- 
ometer, of the blue lobe of the L1157 outflow. We observed four lines at 3 mm, namely 
CH 3 OH (2 k -1k), HC 3 N (11-10), HCN (1-0) and OCS (7-6). Moreover, the bright Bl 
clump has also been observed at better spatial resolution in CS (2-1), CH3OH (2i- 
li)A~, and 34 SO (32~2i). These high spatial resolution observations show a very rich 
structure in all the tracers, revealing a clumpy structure of the gas superimposed to an 
extended emission. In fact, the three clumps detected by previous IRAM-30m single 
dish observations have been resolved into several sub-clumps and new clumps have 
been detected in the outflow. The clumps are associated with the two cavities created 
by two shock episodes driven by the precessing jet. In particular, the clumps nearest 
the protostar are located at the wall of the younger cavity with a clear arch-shape form 
while the farthest clumps have slightly different observational characteristics indicat- 
ing that they are associated to the older shock episode. The emission of the observed 
species peaks in different part of the lobe: the east clumps are brighter in HC3N (11- 
10), HCN (1-0) and CS (2-1) while the west clumps are brighter in CH 3 OH (2 K -1 K ), 
OCS (7-6) and 34 SO (32-2i). This peak displacement in the line emission suggests a 
variation of the physical conditions and/or the chemical composition along the lobe 
of the outflow at small scale, likely related to the shock activity and the precession 
of the outflow. In particular, we observe the decoupling of the silicon monoxide and 
methanol emission, common shock tracers, in the Bl clump located at the apex of the 
bow shock produced by the second shock episode. 

Key words: ISM: individual: L1157 - ISM: jets and outflows - ISM: molecules. 



1 INTRODUCTION 

L1157-mm is a Class protostar located a t 440 pc, with 
a ma ss of ~ 0.12 M© (|Gueth et all Il997l ; iBeltran et ail 
|2004| ). driving a spectacular bipolar outflow. The out- 
flow has been studied in detail through many molecu- 
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With respect to other outflows driven by low-mass pro- 
tostars, L1157 stands out for its rich mm-wave spectrum and 
it can be considered as the prototype of ch emically active 
outflows. IBachiller fc Perez Gutierred (|l997l ) have mapped 
the L1157 outflow in many molecular mm-emission lines 
and six low-velocity molecular clumps have been detected 
along the lobes. These clumps have strong emission due 
to molecular species which increase their a bundances only 
when a shock-induced chemistry is active (Bachiller et al. 
2001). Two main shock events have been detected in the 
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blue lobe of the o utflow and the int erferometric image of 
the CO (1-0) line (|Gueth et al.lll996h reveals the presence 
of two prominent limb-brightened cavities created by the 
propagation of large bow-shocks. The different orientation 
of the two cavities testifies the precession of the driving, 
highly-collimated jet. 

dumpi ness in chemically active outflows has bee n stud- 
ied recently l|Viti et alj|2004 iBenedettini et alj|2006h by fit- 
ting the line profiles of several species by the use of a chem- 
ical and a radiative transfer model. These studies indicate 
that the size derived from single dish observations (~ 0.12 
pc for CB3 and ~ 0.04 pc for LI 157) are upper limits and 
that the clumps contain substructures. In absence of obser- 
vational evi dence of the struc ture of the low velocity c lumps 
in outflow, et al. | <|2004h and IBenedettini et al.1 (|2006l ) 
assumed the density profile derived for starless cores in their 
modelling. However, the outflow clumps show fundamental 
differences with respect to starless cores, for example in the 
temperature (starless cores have temperature ^ 10 K) and in 
the chemical structure (in starless cores most of the molecu- 
lar material is depleted toward the center of the core) . Direct 
evidence of the real structure of these low velocity clumps 
can only be obtained via very high spatial resolution obser- 
vations using an interferometer. 

In order to investigate the small scale structure of the 
low velocity clumps in outflows we carried out interferomet- 
ric observations of the nearby L1157 outflow with the IRAM 
Plateau de Bure Interferometer (PdBI) and the results are 
presented in this paper. We mapped the blue lobe, which is 
the richest in terms of molecular emission, in a sub-sample 
of the species already observed with the IRAM-30m single 
dish in this outflow. In particular, we selected molecules that 
are high density tracers (as CS, HCN and HC3N) and whose 
chemical abundance is significantly affected by the presence 
of a no-dissociative shock (as CH3OH and OCS). The ob- 
servations are described in Sect. 2. We present the results 
and the calculation of the column densities in Sect. 3 and in 
Sect. 4, respectively. The results are discussed in Sect. 5 and 
the main conclusions of the paper are summarized in Sect. 
6. 



2 OBSERVATIONS 

The interferometric observations were carried out with the 
IRAM interferometer at Plateau de Bure in three differ- 
ent periods and configurations. The first set of observations 
was carried out between May and July 1997, the second 
between December 1997 and March 1998 and the third be- 
tween September and November 2004. Some preliminary re- 
sults from the fi r st an d second data sets were presented by 
iPerez-Gutierrezl <|l999T ). but we report here final images to- 
gether with a detailed analysis of the three data sets. In 
the first two periods the CD configuration was used with 
three configurations of 4 or 5 antennas (5C2+4D1+5D) and 
a baseline between 24 m and 176 m, with a final clean beam 
of 6.7x5.5 arcsec 2 at 85 GHz and 3.8x3.5 arcsec 2 at 100 
GHz. A mosaic of 5 fields was performed to cover all the 
south lobe of the L1157 outflow at 3 mm. The correlator 
units were configured in order to have a resolution of 78 kHz, 
which corresponds to 0.25 km s _1 at 3 mm. In a first setup 
we observed CH3OH (2 K -1 K ) at 96.741 GHz in the LSB 



and HC 3 N (11-10) at 100.076 GHz in the USB and in a sec- 
ond observation HCN (1-0) at 88.632 GHz and OCS (7-6) 
at 85.139 GHz were observed simultaneously. The CH3OH 
(5k-4k) line at 241.285 GHz was observed simultaneously. 
The phase and amplitude calibration was achieved by ob- 
servations of 2021+614, which is close to L1157 in the sky. 
The band pass of the receivers were calibrated by observa- 
tions of 3C273 and 3C454.3. To correct for decorrelation 
due to phase noise, amplitude calibration was also done 
relative to 2021+614, whose flux was determined relative 
to 3C273 and 3C454.3. In the observations carried out in 
2004 the CD configuration was used with 6 or 5 antennas 
(5D+6Cp) obtaining a better spatial resolution of ~ 3.1x2.8 
arcsec 2 at 97 GHz. One pointing on the central zone of the 
blue lobe, the Bl clump [R.A.(J2000) 20 h 39 m 9 ! :5; dec(J2000) 
68°01'15"], was performed to observe simultaneously CS (2- 
1) at 97.981 GHz, CS (5-4) at 244.935 GHz and CH3OH 
(5 K -4 K ) at 241.791 GHz. The 34 SO (3 2 -2i) at 97.715 GHz 
and CH3OH (2i-li)A" at 97.583 GHz were automatically 
covered by the setting of the autocorrelator units that were 
configured in order to have a resolution of 0.5 km s _1 at 
3 mm. The phase and amplitude calibration was achieved 
by observations of 2037+511 and 1928+738, which are close 
to L1157 in the sky. The band pass of the receivers was 
calibrated by observations of 2145+067 and 1928+728. The 
flux calibration was determined relative to MWC349. The 
data were calibrated and analyzed with the GILDAS soft- 
ware, the images were produced using natural weighting and 
cleaned in the usual way. The details of the observations are 
summarized in Table [T] The 1.3 mm data (not shown) have 
a very low signal to noise ratio (S/N), this does not allow 
us to properly analyse those data. 

In the interferometers the emission from structures 
more extended than 2/3 of the primary beam is filtered 
out. The extended emission, that is expected to be present 
in outflows, could be partially recovered by adding single 
dish data to the interferometric observations. To this aim, 
we include the short-spacings to the PdBI data for the two 
lines CH3OH (2 k -1k) and HC N (1-0), for which sin gle dish 
IRAM-30m data are available (|Bachiller et alj|200ll ). 



3 RESULTS 

3.1 Overall structure 

The maps of the blue lobe of the L1157 outflow in HCN 
(1-0), HC 3 N (11-10), CH3OH (2 k -1k) and OCS (7-6) are 
shown in Fig.Q] The central part of the outflow has also been 
mapped in three more lines, namely CS (2-1), CH3OH (2i- 
li)A" and 34 SO (3 2 -2i) whose maps are shown in Fig. [2] 
These high spatial resolution observations reveal the pres- 
ence of a clumpy structure superimposed to an extended 
emission. 

The extended emission is clearly visible in the two 
maps, HCN (1-0) and CH3OH (2 k -1k), where the short- 
spacings are included, while it is not present in the other 
maps, where it is probably filtered out by the interferom- 
eter. The extended emission matches very well with pre- 
vious observations, n amely multi-species single dish maps 
llBachiller et al.l l200ll ) and interferometric CO (1-0) map 
IjGueth et al.lll996P l. tracing the large-scale outflow struc- 
ture. 
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Table 1. List of the observed transitions and observing parameters. The first three lines were observed only in the central part of the 
blue lobe, the Bl clump. 



transition 


V 


resolution 


PdBI 




clean beam 


rms 


short-spacings 




GHz 


km s — 1 


configuration 




arcsec 


Jy beam -1 km s -1 




CH3OH (2i-li)A" 


97.583 


0.5 


5D+6Cp 


3.14 


X 2.82 (P.A.=69°) 


0.003 


N 


34 SO (3 2 -2i) 


97.715 


0.5 


5D+6Cp 


3.14 


x 2.81 (P.A.=69°) 


0.003 


N 


CS (2-1) 


97.981 


0.5 


5D+6Cp 


3.11 


x 2.79 (P.A.=70°) 


0.009 


N 


OCS (7-6) 


85.139 


0.28 


5C2+4D1+5D 


6.68 


x 5.50 (P.A.=-79°) 


0.05 


N 


HCN (1-0) 


88.632 


0.27 


5C2+4D1+5D 


5.63 


x 4.28 (P.A.=93°) 


0.80 


Y 


CH3OH (2 K -1 K ) 


96.74lt 


0.24 


5C2+4D1+5D 


5.06 


x 4.37 (P.A.=-85°) 


1.30 


Y 


HC 3 N (11-10) 


100.076 


0.23 


5C2+4D1+5D 


3.84 


x 3.51 (P.A.=-43°) 


0.14 


N 
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Figure 1. Intcrfcromctric images of the blue lobe of the L1157 outflow. Note that for HCN (1-0) and CH3OH (2k-1k) the single dish 
data have been added to the interferometric data. The star indicates the position of the millimeter driving source and the ellipse at the 
bottom right shows the clean beam. The contours are: for HCN (1-0) first contour 2.5 Jy beam -1 km s -1 (3a), level steps 2 Jy beam -1 
km s -1 ; for CH3OH (2k-1r) nrs t contour 4 Jy beam -1 km s -1 (3a), level steps 1.5 Jy beam -1 km s -1 ; for HC3N (11-10) first contour 
0.4 Jy beam -1 km s -1 (3a), level steps 0.15 Jy beam -1 km s -1 and for OCS (7-6) first contour 0.15 Jy beam -1 km s -1 (3a), level 
steps 0.05 Jy beam -1 km s -1 . 




20 h 39 m 10 s 
R.A. 

Figure 2. Maps of the Bl clump of the L1157 outflow. The contours are: for CS (2-1), first level is 0.027 Jy beam -1 km s -1 (3(t), level 
steps are 0.02 Jy beam -1 km s -1 ; for CH3OH (2i-li)A - , first level is 0.009 Jy beam -1 km s -1 (3ct), level steps are 0.005 Jy beam -1 
km s -1 ; for 34 SO (3 2 -2i), first level is 0.01 Jy beam -1 km s -1 (3<r), level steps are 0.004 Jy beam 1 km s 1 . The ellipse at the top 
right shows the clean beam. The 34 SO map has been convolved at a higher spatial resolution with respect to the observations, in order 
to improve the S/N. 
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Table 2. Coordinates of the clumps. In the last column the tran- 
sition used to define the coordinates is listed. 



clump R.A.(J2000) Dcc(J2000) transition 





h 


m 


s 




O 


/ 


II 




BOa 


20 


39 


10. 


.4 


68 


01 


38 


HCN (1-0) 


BOb 


20 


39 


05. 


.7 


(58 


01 


12 


CH 3 OH (2 K -1 K ) 


BOc 


20 


39 


09. 





(58 


01 


3(5 


CH3OH (2 K -1 K ) 


BOd 


20 


39 


0V. 


7 


(58 


01 


28 


CH3OH (2 K -1 K ) 


BOe 


20 


39 


10. 


3 


68 


01 


24 


CH3OH (2 K -1 K ) 


Bla 


20 


39 


10. 


2 


(58 


01 


12 


HCN (1-0) 


Bib 


20 


39 


08. 


.8 


(58 


01 


10 


CH3OH (2 K -1 K ) 


Blc 


20 


39 


10. 


.4 


(58 


01 


07 


CS (2-1) 


Bid 


20 


39 


09. 


2 


(58 


01 


04 


34 SO (3 2 -2i) 


B2c 


20 


39 


11. 


.8 


68 


00 


58 


CH3OH (2 K -1 K ) 


B2a 


20 


39 


12. 


6 


(58 


00 


10 


CH3OH (2 K -1 K ) 


B2b 


20 


39 


14. 


1 


68 


00 


24 
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The interferometric observations show also a very rich 
structure at small-scale, revealing a dumpiness of the molec- 
ular gas richer than the one previously observed. The three 
large molecular clumps, called B0, B l and B2, previously 
detected by single dish observations l|Bachiller et alj|200ll ) 
have been resolved into more sub-clumps: the B0 clump is re- 
solved into three sub-clumps labelled BOa, BOc and BOe, the 
Bl clump is resolved into four sub-clumps labelled Bla, Bib, 
Blc and Bid, and the B2 clump is resolved into two sub- 
clumps labelled B2a and B2b. Moreover, three new clumps 
have been detected: BOb at north-west, BOd located between 
BOb and Bl and B2c located between Bl and B2. The co- 
ordinates of the identified clumps are listed in Table [2] they 
have been defined by using the tracer where the clump is best 
defined. Some of these clumps are associated with the shock 
front where the outflow jet impacts the ambient medium 
while others are in the wake of the shocked region where the 
matter has already been affected by the passage of the shock 
(see Sect. 15.21 for a detailed description). 

In Fig. |3] we show the spectra of HCN (1-0), CH3OH 
(2 k -1k), HC 3 N (11-10) and CS (2-1) in some of the iden- 
tified clumps. At the spectral resolution of the observations 
the three lines of the hyperfine structure of the HCN (1-0) 
transition at 88.630 GHz, 88.632 GHz and 88.634 GHz are 
resolved. The three lines are strongly self-absorbed along the 
flow but not at the position of the mm source. The velocity of 
the self-absorption features remains constant along the lobe 
suggesting the presence of an extended ambient component 
of HCN; alternatively the self-absorption might be due to 
high opacity effects. In the spectra of the CH3OH (2k-1k) 
transition (see Fig. [3J three lines are detected: (2_i-l_i)E 
at 96.739 GHz, (2 -l )A+ at 96.741 GHz and (2 -l )E at 
96.745 GHz. The three lines are blended in most of the map. 

The clumps are not detected in all the species. In fact, 
the S/N of the maps is not the same due to the fainter signal 
of some lines. In particular, the maximum S/N in the map 
is higher for HCN (1-0), CH3OH (2 k -1k) and CS (2-1) (19, 
18 and 16 respectively) while it is lower for HC3N (11-10), 
OCS (7-6) and 34 SO (3 a -2i) (7, 5 and 4 respectively). 

The CH3OH (2 k -1k) and HCN (1-0) emission shows 
the presence of an extended molecular component (note that 
these maps include the short-spacings) and they also have 
the richest structure in term of clumps (see Fig. [1} • In par- 
ticular, all the identified clumps, but BOb and BOc, are de- 



tected in these tracers. BOb and BOc are detected only in the 
methanol lines. The HC3N (11-10) emission has a definitely 
lower S/N; nevertheless 5 clumps are clearly identified (see 
Fig. [T]). In the northern region of the map the peak of the 
HC3N (11-10) emission is at an intermediate position be- 
tween the BOa and BOe clumps identified by HCN. The Bla 
clump is brighter with respect to the Bib. The two clumps 
B2a and B2b have also been detected in HC3N; however, 
the B2b peak is shifted toward the north-west with respect 
to the peak of HCN (1-0) and CH3OH (2 k -1k). Despite 
the low S/N of the OCS map (see Fig. [TJ, the two most 
southern clumps B2a and B2b are clearly identified, while 
only a weak signal can be seen in BOd and Bib. The quality 
of the 34 SO map (see Fig. [2)| is also poor, however a clear 
arch-shape emission is detected. Three sub-clumps are de- 
tected above the noise: Bla and Bib have also been seen 
in other tracers, while a new clump, Bid is detected at an 
intermediate position between Bib and Blc. This clump is 
clearly defined only in this tracer, however, given the low 
S/N of the peak, we cannot rule out the possibility that it is 
an artifact. Despite the low S/N, it seems that hints of 34 SO 
and OCS emission are mainly peaking towards the western 
part of the Bl clump and in the southern part of the blue 
lobe. 

Finally, note that most of the clumps are quasi-circular 
with size of 10 - 15 arcsec (0.02 - 0.03 pc at a distance of 440 
pc), apart from B2a and B2b that are elongated with size of 
~ 10x20 arcsec 2 . The emission in the most southern clumps 
is at higher velocity with respect to the northern clumps. 
This trend can be seen in the channel map of HCN in Fig. 2] 
and it is also present in all the other species (channel maps 
not shown). 

3.2 The clump BO 

BOa, BOb, BOc, BOd and BOe are the clumps nearest to the 
protostar driving the outflow. BOa can be identified with the 
peak of the B0 clump observed in the single dish observa- 
tions. Close to BOa we detect other clumps: BOc and BOd at 
south-west of BOa, BOe at south, and BOb at west. BOb and 
BOc are detected only in CH3OH while BOa is detected also 
in the HCN (1-0) (note that these clumps are not mapped 
in CS and 34 SO). BOd is detected in HCN (1-0) and CH 3 OH 
and emits marginally in OCS (7-6) and CS (2-1), while BOe 
is detected in all the tracers but OCS and 34 SO. The HC 3 N 
(11-10) line shows an elongated emission covering BOa and 
BOe with a peak at an intermediate position between the 
two. 



3.3 The clump Bl 

Bl is the brightest clump of the outflow and it has been 
mapped in CH 3 OH (2i-li)A~, CS (2-1) and 34 SO (3 a -2i) 
with a spatial resolution higher than the other maps. Bl has 
a well defined arch-like shape with four sub-clumps (Bla, 
Blc, Bid and Bib) outlining the arch. The CH3OH (2i- 
li)A - emission shows the same quasi-circular morphology 
of the larger CH3OH (2r-1k) map with a peak in Bib while 
the CS emission shows three clumps located in an arched- 
like shape with Blc showing the strongest emission. The 
arch shape emission is also detected in the HC3N (11-10) 
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Figure 3. Spectra of the HCN (1-0), CH 3 OH ((2 K -1 K ), HC3N (11-10) and CS (2-1) in the BOa, BOd, BOe, Bla, Blc, Bib, and Bid 
clumps. The Tmb of the HC3N line is multiplied by a factor of 5. The dashed line is the t>i sr =2.7 km s . The continuous lines mark 
the position of the three hyperfine transitions of HCN (1-0) and the three transitions CH3OH (2_j— l_i)E, (2q— lrj)A+ and (2rj— lo)E. 



and 34 SO (3 2 -2i) lines. In CS (2-1) we detect a diffuse com- 
ponent that is elongated in the so uth direction , simila rly to 
what is seen in the SiO map by iGueth et al, I i|l998l ). The 
elongated feature is clearly visible in the channels map (Fig. 
for velocities higher than —5 km s . The channel map 
also shows that close to the systemic velocity («i 8r = 2.7 
km s -1 ), the diffuse emission assumes a more pronounced 
arch-like shape and that the BOe clump becomes elongated 
toward Bla and then it 'connects' to the Bl clump at v = 
— 1.4 km s _1 . In the Bla position a second CS (2-1) veloc- 
ity component is identified in addition to the main compo- 
nent centered at km s _1 (see Fig. [3}; the second component 
peaks at —11. 4 km s _1 and extends up to —16 km s _1 . This 
component is a l so pre sent in the SiO (3-2) line observed by 
iBachiller etafl (|200lh . In the channel map of the CH3OH 
(2i-li)A~ line (shown in Fig. [6j one can see that the peak 
position of the Bib clump is shifted toward the south at 
higher velocity. Indeed, the emission peaks at an offset of 
(— 4",4") with respect to the center of the map in the —5.4 
km s _1 channel and at an offset (—4",— 2") in the 2.6 km 
s _1 channel. 

3.4 The clump B2 

In the southern zone of the blue lobe (not mapped in CS and 
,i4 SO), the interferometric observations show the presence of 
three clumps. The faintest clump is B2c, detected only in 
HCN (1-0) and CH3OH (2 K -1 K ). It is located between Bl 
and B2 and its shape and size are more similar to the north- 
ern clumps than the other two B2 sub-clumps. In fact, B2a 



and B2b show observational characteristics slightly different 
with respect to the other clumps. They are detected in all the 
molecules, with very similar intensity. In all transitions B2a 
and B2b have a line peak velocity (v ~ 2 km s" 1 ) closer to 
the ambient velocity vi SI = 2.7 km s -1 while all other clumps 
have a typical peak velocity of ~ km s~ . The final veloc- 
ity of the blue line wing is smaller of about 2-4 km s _1 in 
B2a and B2b than in the other clumps. B2a and B2b are the 
clumps farthest from the protostar and also their morphol- 
ogy is different: they are not quasi-circular or displaced in an 
arch-like shape but they have an elongated shape of about 
10 x 20 arcsec 2 and are displaced along the outflow axis. All 
these characteristics suggest that B2a and B2b are associ- 
ated with an older shock episode and they have been slowed 
down during their propagation in the molecular cloud. The 
lower velocity of the B2 clumps might also be caused by pro- 
jection effect related to the precession of the outflow axis. 



4 COLUMN DENSITIES 

We derived the molecular column density in the detected 
clumps, from the integrated intensity of the observed emis- 
sion lines, assuming LTE condition and that the lines are 
optically thin. In this case the following formula can be used 

N = 1.67 x 10 14 exp ( [ Tmb dv (1) 

where Q(T Iot ) is the partition function, \x is the dipole mo- 
ment in debyes, v is the frequency in GHz, S is the line 
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Figure 4. Channel map of the blue lobe of the L1157 outflow in the HCN (1-0) line at 88.632 GHz. The thick box points out the 
ambient velocity emission. First level and level steps are 0.2 Jy beam -1 (3<r). 



strength, E up is the energy of the upper level of the tran- 
sition and the integral of the line emission is in K km s" 1 . 
The integration limits used to calculate the column density 
is the same for all molecules in each clump and it has been 
determined using the molecule where the clump is best de- 
fined. The column densities for the observed molecules are 
reported in Table [3] These values must be considered as 
lower limits because: i) the LTE assumption may be not 
valid for most of the transitions since they have a critical 
density ^10 5 cm" 3 ; and ii) the assumed temperature of 60 
K for B 2a and B2b and 80 K fo r the others clumps, de- 
rived bv lTafalla fc Bachillerl (|l995l ) by means of VLA obser- 
vations of NH3 emission, are lower limits since at these high 
temperatures the am monia rotational tra nsitions are prob- 
ably not thermalized (|Danbv et al.|[l988l ) and the real gas 
kinetic temperature may be significantly higher. Moreover, 
the HCN column density is significantly underestimated be- 
cause of the strong self-absorption of the observed line. For 
the calculation of the CH3OH column density we used the 
(2i-li)A~ line at 97.583 GHz in the clumps where it has 
been observed (i.e. BOd, BOe, and the Bl sub-clumps). For 
the other clumps the situation is more difficult because we 
must use the (2k-1k) data and the three lines at 96.739, 
96.741 and 96.745 GHz are closely blended (see Fig.©. How- 



ever, to roughly evaluate the methanol column density in 
the whole lobe we used the most intense line (2o~lo) A + at 
96.741 GHz. The integrated line intensity for this line has 
been calculated summing the emission in a velocity range 
that has been defined by eye in each clump in order to min- 
imize the contribution of the adjacent lines. However, we 
stress that in most cases this contribution is not negligible, 
leading to an overestimate of the methanol column density 
as it can be seen by comparing the values derived using this 
line with those derived by using the (2i-li)A~ line in the 
clumps where it is observed (column 5 vs column 6 of Ta- 
ble [3]l : the first values are always higher (never more than a 
factor of 10) than the second. 

The CS and HC3N column densities are quite uniform 
in all the clumps, (3-6) xlO 14 cm" 2 and (1-3) xlO 13 cm" 2 , 
respectively. Note that the CS column density in the Bla 
clump is the sum of the two velocity components detected 
in this position. On the other hand, the methanol column 
density shows a higher value, 1.2 xlO 16 cm -2 , in Bib with 
respect to the typical value of the other clumps, (2-5) x 10 15 
cm" 2 . On the contrary, the HCN column density has one of 
the lowest values 5xl0 14 cm -2 in Bib while it is higher in 
Bla and Blc, ~ 2xl0 15 cm" 2 . The 34 SO column density is 
~ 3xl0 13 cm" 2 in Bla and Blc while it is slightly higher, 
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(8-9) xlO 14 cm~ 2 , in Bib and Bid. In all the species, the 
column densities in B2a and B2b have similar values and are 
lower than the values in Bl. Note, however, that the size of 
the clumps where the column densities are calculated are 
larger in B2a and B2b than in the other clumps. 

The column densities derived by merging our PdBI data 
with previously obtained IRAM-30m data for HCN and 
CH3OH are similar (difference of 10% on average) to the 



values derived bv lBachiller fc Perez Gutierrez! l| 19971 ) in Bl 
and B2 using only single dish IRAM-30m data. This shows 
that most of the large scale emission filtered out by the in- 
terferometer has been recovered adding the single dish data 
to the interferometric data. 
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Table 3. Molecular column densities in the clumps. 



clump 


N(HC3N) 
(cm -2 ) 


N(HCN) 
(cm -2 ) 


N(OCS) 
(cm -2 ) 


N(CH 3 OH)t 
(cm -2 ) 


N(CH 3 OH)t 
(cm ) 


N( 34 SO) 
(cm -2 ) 


N(CS) 
(cm -2 ) 


BOa 


2.0xl0 13 


1.3xl0 15 




9.7xl0 15 








BOb 








1.7xl0 15 








BOc 








9.7xl0 15 








BOd 




2.4xl0 13 


6.4xl0 13 


1.2xl0 16 


3.6xl0 15 




9.3xl0 13 


BOe 


1.7xl0 13 


l.lxlO 15 




5.5xl0 15 


2.3xl0 15 




5.3xl0 14 


Bla 


3.1xl0 13 


1.8xl0 15 




3.5xl0 16 


3.4xl0 15 


3.6xl0 13 


5.5xl0 14 


Bib 


3.4xl0 13 


5.0xl0 14 


1.5xl0 14 


5.0xl0 16 


1.2xlO ie 


8.3xl0 13 


4.1xl0 14 


Blc 


3.4xl0 13 


1.7xl0 15 




4.0xl0 16 


3.5xl0 15 


2.6xl0 13 


6.4xl0 14 


Bid 


2.2xl0 13 


9.3xl0 14 




3.8xl0 16 


5.0xl0 15 


9.0xl0 13 


3.1xl0 14 


B2c 




6.3xl0 13 




9.3xl0 15 








B2a 


1.7xl0 13 


4.7xl0 14 


7.9xl0 13 


1.7xl0 16 








B2b 


9.2xl0 12 


4.4xl0 14 


7.1xl0 13 


1.7xl0 16 








t derived from the (2o-lo)A+ line 
t derived from the (2i-li)A — line 



5 DISCUSSION 

5.1 The dumpiness of the blue lobe 

The two species for which the single dish data have 
been added to the interferometric data, namely HCN and 
CH3OH, show the presence of a diffuse gas component that 
matches the morphology of the blue lob e observed in other 
species such as CO CS, SiO and SO (|Gueth et al J 1 19961 ; 
iBachiller et~aH200ll i. Our high spatial resolution maps show 
that the outflow emission is indeed continuous, with a num- 
ber of local peaks. In fact, the large clumps detected by pre- 
vious single dish observations have been resolved into sev- 
eral smaller clumps and few new clumps have been detected 
along the lobe. T hese findings confirm what suggested by 
previous studies l|Viti et al. 1 12004| ; iBenedettini et all 12009 ) 
that several gas components are present in the low veloc- 
ity molecular clumps observed in chemically active outflows 
with single dish telescopes. 

Two main structures are identified (see Fig. [T}: the 
southern one is associated with the B2 clumps and the north- 
ern one is associated with the BO and Bl clumps. These two 
structures of diffuse molecular gas are associa ted with the 
two ca vities identified in the CO (1-0) map of lGueth et al.l 
(1996). The two cavities have probably been created by the 
propagation of large bow-shocks due to a highly collimated, 
precessing jet. It is worth noting the two continuous struc- 
tures perpendicular to the flow axis: the first connecting 
BOa, BOc and BOd and the second slightly to the north at 
the level of BOb; they seem to be detected in HCN and 
CH3OH as well as CO. The nature of these structures re- 
mains unclear. 

In Fig.[7]we show the CS (2-1) and CH3OH (2i-li)A - 
maps superimposed to the CO (1-0) map. The clumps iden- 
tified in our maps are located at the walls of the C2 cavity. 
In particular, BOa and BOe trace the east wall, BOb and 
BOd trace the west wall and the arch-shaped Bla, Bib, Blc 
and Bid clumps are located at the apex of the cavity. Our 
maps support the scenario that the cavity was created by 
the propagation of large bow-shock, showing the presence of 
clumpy structures at the walls of the cavity with the arch- 
shaped clumps at the apex that may trace the front of the 
bow-shock. The situation in the older CI cavity is less clear 
since this cavity is less defined in the CO map. However, it 



seems that the B2a and B2b clumps trace the west wall of 
the cavity. 



5.2 East /West chemical difference 

The observed species peak in different parts of the C2 cav- 
ity. Considering also other interferometric maps found in 
the literature [SiO (2-1), spatial resolution ~ 2.5 arcsec 
|Gueth et al.lll998l) and NH 3 (3,3), ( 1,1), spatial resolution 
~ 5 arcsec ~(|Tafalla fc Bachillerlll995l )]. we can identify two 
groups of molecules, one peaking in the west-side and the 
other peaking in the east-side. In particular, the east clumps 
BOa, BOe, Bla and Blc are brighter in HC 3 N (11-10), HCN 
(1-0), CS (2-1) and NH 3 (3,3) and SiO (2-1) while the west 
clumps BOb, BOd and Bib are brighter in CH3OH (2 k -1k), 
OCS (7-6) and 34 SO (3 2 -2i). This peak displacement in the 
line emission suggests an inhomogeneity in the physical and 
chemical conditions of the observed clumps. 

The different behaviour between the east and west parts 
of the map is also clear from the column densities (see 
Table |3J. N(CS) and N(HC 3 N) are quite constant in all 
the clumps; N(HCN) is higher in the east clumps while 
N(CH 3 OH), N(OCS) and N( 34 SO) are higher in the west 
clumps. In particular, the OCS and CH3OH column density 
shows the highest value in Bib while here the HCN column 
density assumes its lowest value. The 34 SO column density 
is slightly higher in Bib and Bid than in Bla and Blc. 
Despite the large uncertainties associated with the derived 
column densities, these differences suggest that a chemical 
differentiation is present among the clumps associated with 
the C2 cavity. In particular, the N(CH 3 OH)/N(CS) ratio 
is higher in the west part (38 in BOd, 29 in Bib and 16 
in Bid) than in the east part (6 in Bla and Blc and 4 in 
BOe) (see Fig.0. Similarly N( 34 SO)/N(CS) is higher in the 
west (0.20 in Bib and 0.29 in Bid) than in the east (0.06 in 
Bla and 0.04 in Blc). Since the CS column density is quite 
constant throughout the map, the difference in the ratio is 
an indication of an higher chemical abundance of methanol 
and sulfur oxide in the west part of C2 cavity, especially 
in Bib for CH3OH and in Bid for 34 SO. This result is in- 
dependent on the assumption that the temperature of 80 
K is constant throughout the Bl clump. In fact, 80 K was 
derived from the non uniform NH3 (3,3)/(l,l) line ratio by 
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Figure 7. Left panel: overlay of CO (1-0) (grey scale) jGueth et al.lll996l) , CS (2-1) (black line) and CH 3 OH (2i-li)A - (white line). 
The star marks the position of the protostar and the two straight lines indicate the direction of the precessing jet. Right panel: ratio of 
the CS and CH3OH (2i-li)A" emission; contours are: 10, 20, 30, 40, 50, 60, 80, 100, 150, 200, 250, 300. 



iTafalla fc Bachillerl l|l995l ). In particular, the reference value 
of 80 K is suitable for the east clumps where the ammonia 
peaks but it may not be suitable for methanol, which peaks 
in the west. However, by changing the temperature from 
10 K to 200 K, the methanol column density changes only 
by a factor of 3, and the N(CH 3 OH)/N(CS) ratio is always 
significantly higher at west than at east. 



5.3 The Bl internal structure 

The asymmetry in the line emission is particularly evident in 
Bl. In Fig|8]we compare the emission of CS (2-1), CH3OH 
(2i-li)A~ and SiO (2-1 ). CS and SiO, as well as NH 3 
jTafaUa fc Bachillerl Il995h , show a similar behaviour, peak- 
ing ~ 10 arcsec east of the peak of CH3OH. Since the Bl 
clump is located at the apex of the bow shock produced 
by the second shock episode, it is particularly interesting 
to note the different spatial distribution of the SiO, NH 3 , 
and CH3OH molecules that are common shock tracers. SiO 
(2-1) and CH3OH (2k-1k) are expected to be excited by 
the same physical conditions since their critical density and 
excitation temperature are similar (n CT it ~ 10 6 cm* 3 and 
T cx ~ 4-15 K). It is then probable that a different chem- 
istry is at the origin of the observed displacement, although 
small differences of the gas density and temperature are also 



possible between the subclumps. In fact, me thanol and am- 
moni a are abundant on grain ice mantles (|Charnlev et al.l 
1992) and they evaporate in th e gas phase as the t emper- 
ature increases above ~ 120 K (|Collings et alj|20040 . while 
SiO is rapidly formed in the gas phase after the sputter- 
ing of significant amounts of silicon from th e dust grains 
in C-shock with speeds larger than 25 km s -1 l|Schilke et alj 
1 19971 ; IP ineau des Forets et al.|[T997h . Since the efficiency of 
the Si sputtering increases with the velocity of the shock, 
the higher brightness of the SiO emission in the eastern 
part of the C2 cavity might suggest that the velocity of 
the shock is higher in the east wall of the cavity with re- 
spect to the west wall. This is also supported by the fact 
that, in general, the terminal velocity of the blue wing for 
the observed lines is higher, from 4 up to 10 km s~ , for the 
eastern clumps than for the western clumps (see Fig. [3} . On 
the other hand, the chemical models developed for th e low- 
velocity molecular clumps in outflows l|Viti et al.ll2004l ) show 
that the methanol chemical abundance is higher in models 
with high gas density. In particular, the abundance increases 
by more than one order of magnitude as the gas density in- 
crea ses from 10 5 cm " 3 to 10 6 cm -3 (see Table 5 and Fig. 
3 of lyitT et al.l |2004| ) . This suggests that the west clumps 
might have a slightly higher density, nH 2 ~ 10 6 cm* 3 , than 
the mean value hh, = 3xl0 5 cm* 3 derived in Bl from the 
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Figure 8. Overlay of the SiO (2-1) (grey scale) jGueth et alj|l998ft with: a) CH 3 OH (2i-li)A - (continuous line) and b) CS(2-1) 
(continuous line). For SiO the first level is 1 Jy beam -1 km s — 1 (3<r), level steps are 0.5 Jy beam -1 km s — 1 ; for CH3OH (2i— li)A - , first 
level is 0.01 Jy beam -1 km s -1 (3cr), level steps are 0.007 Jy beam -1 km s — 1 ; for CS (2-1), first level is 0.027 Jy beam -1 km s — 1 (3<r), 
level steps are 0.02 Jy beam -1 km s -1 . 



LVG fitting of the SiO lines observed w ith the IRAM-30m 
and JCMT telescopes (jNisini et alj|2007h . 

The asymmetry in the emission of the shock tracing 
species may be explained in the framework of the morphol- 
ogy of the blue lobe of the L1157 outflow. The lobe is com- 
posed by two big cavities created by two subsequent bow 
shocks driven by a precessing jet (|Gueth et al.lll996l ): the 
first shock episode is associated with the CI cavity and the 
southern clumps and the second shock episode is associated 
with the C2 c avity and the northern clumps. The model by 
iGueth et al.l (|l996l ) found that the jet precesses on a cone 
with opening angle of 6° and that the CI cavity lies mainly 
on the plane of the sky while C2 is elongated towards us 
and towards the west. They also found that between the 
two shocks the jet precessed a quarter of the rotation pe- 
riod. Due to the precession of the jet axis from east to west, 
the west wall of the C2 cavity is expanding into the denser 
medium of the harboring molecular cloud while the east wall 
is expanding in a medium already swept-up by the older 
shock episode. Therefore chemical differentiations between 
the east and west clumps are expected because of the differ- 
ent physical conditions of the pre-shock gas. In this scenario 
one would expect a higher SiO abundance at east, where we 
observe a higher velocity, and higher CH3OH abundance at 
west, where the density is likely to be higher. On the other 
hand, the change of the outflow inclination with respect to 
the plane of the sky between the two shock episodes, even 
if quite small (~ 6°), could imply that the eastern part of 
Bl may intersect the dense wall of the previously formed 
CI cavity. If this is the case, one would expect a higher 



gas density at the east of Bl, i.e. an opposite scenario with 
respect to the previous one. A similar puzzling finding is 
observable in the south part of Bl where a linear structure 
of diffuse emission is detected in CS (2-1) and in SiO (2-1) 
but not in CH3OH (2i-li)A~ (see Fig. E]). Also in this case 
CS and SiO trace the same gas comp onent and they are 
decoupled from the methanol emission. IGueth et al.l (|l998l ) 
suggested that the SiO finger, that lies along the jet direc- 
tion, traces the magnetic precursor of the shock. An en- 
hancement of the SiO abundance due to the interaction of 
the magnetic and/or radiative precursor of the shock with 
the amb ient gas has also been ob served in the L1148-mm 
outflow (|jimenez-Serra et al.ll2004h . 



Finally, it is worth noting that a displacement between 
the CH3OH (2i-li) and SiO (2-1) emission has been also 
observed in the out flow driven by the h igh-mass protostar 
IRAS 20126+4104 (|Cesaroni et al.ll2005l . see their Fig. 16). 
Also in this case, the authors interpret such displacement by 
the light of the outflow precession with the methanol abun- 
dance increased as a consequence of the entrainment of the 
surrounding high density cloud. In conclusion, the real na- 
ture of the observed displacement of the shock tracer species 
in Bl can be explained only with a proper determination of 
the actual physical conditions (density as well as tempera- 
ture) of the clumps. In particular, to this aim observations 
of higher excitation CS transitions are needed to constrain 
possible density variations. 
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6 CONCLUSIONS 

High spatial resolution maps of the blue lobe of the L1157 
outflow, obtained with the PdBI have been presented. The 
blue lobe has been mapped in CH3OH (2 k -1k), HC3N (11- 
10), HCN (1-0) and OCS (7-6) at a spatial resolution of ~ 
6x5 arcsec 2 . Moreover, the bright Bl clump has also been 
observed in CS (2-1), CH3OH (2i-li)A _ , and 34 SO (3 2 -2i) 
at higher spatial resolution of ~ 3x3 arcsec 2 . 

These observations show a very rich structure in all the 
tracers, revealing a clumpy structure of the gas, richer than 
the one observed in previous single dish observations. We de- 
tect new clumps, labeled BOb, BOd, and B2c. Moreover, the 
previously known BO, Bl and B2 clumps have been resolved 
into three (BOa, BOc, BOe), four (Bla, Bib, Blc and Bid) 
and two (B2a and B2b) sub-clumps, respectively. These re- 
sults confirm that the shocked molecular gas in chemically 
active outflows is mainly structured in small clumps with 
size of the order of 0.02 - 0.04 pc. 

The two southern clumps (B2a and B2b) of the blue 
lobe of the L1157 outflow appear similar both in the line in- 
tensity and the elongated shape (20x10 arcsec 2 ). They show 
observational properties slightly different with respect to the 
other clumps, probably due to their older age. The northern 
clumps (BOa, BOb, BOc, BOd, BOe, Bla, Bib, Blc and Bid) 
have typical size of 10 - 15 arcsec and they are located at 
the walls of the C2 cavity formed by the propagation of a 
bow shock, with the Bl sub-clumps tracing the apex of the 
bow shock and the B0 sub-clumps located in the wake of the 
shocked region where the matter has already been affected 
by the passage of the shock. The observed species peak in dif- 
ferent parts of the cavity. By taking into consideration other 
interferometric maps found in the literature, we identify two 
groups of molecules: HC 3 N, HCN, CS, NH 3 and SiO peak 
at the west-side and CH3OH, OCS and 34 SO peak at the 
east-side of the cavity. This peak displacement is also clear 
from the column densities and suggests an inhomogeneity in 
the physical conditions and/or chemical composition of the 
observed clumps. The observed asymmetry may be related 
to the precession of the outflow: because of the precession 
of the outflow axis from east toward west, the west wall of 
the C2 cavity is expanding into the denser medium of the 
harboring molecular cloud so the velocity of the shock may 
be lower and the density of the clumps higher than at the 
east wall which is expanding in a medium already swept-up 
by the older shock episode. However, other hypotheses can- 
not be ruled out, therefore the determination of the nature 
and the physical conditions of the clumps require further 
high spatial resolution observations and the use of a detailed 
chemical and shock model. 

Finally, we note that the shocked gas in the L1157 out- 
flow will be extensively studied by the PACS and HIFI 
instruments on board the Herschel satellite, that will be 
launched in 2008. In particular, the blue lobe is the cho- 
sen position of an unbiased spectral survey in the ~ 450- 
1990 GHz window with a spatial resolution ranging from ~ 
9 to 46 arcsec. Even if the Herschel observations will be at 
lower spatial resolution with respect to the interferometric 
data presented here, the multiline analysis of data from the 
millimeter to the far-infrared spectral range will allow to 
trace all the gas excitation from few tens to few thousands 
of Kelvin and to fully constraint the physical properties of 



the warm gas components along the chemically active L1157 
outflow. 



ACKNOWLEDGMENTS 

We acknowledge Riccardo Cesaroni and Paola Caselli for 
helpful discussions. This work is based on observations car- 
ried out with the IRAM Plateau de Bure Interferometer. 
IRAM is supported by INSU/CNRS (France), MPG (Ger- 
many) and IGN (Spain). For the data reduction MB has 
benefited from research funding from the European Com- 
munity's Sixth Framework Programme. SV acknowledges fi- 
nancial support from an individual PPARC Advanced Fel- 
lowship. 



REFERENCES 

Avery L.W., Chiao M., 1996, ApJ, 312, 788 
Bachiller R., Perez Gutierrez M., 1997, ApJ, 487, L93 
Bachiller R., Perez Gutierrez M., Kumar M.S.N. , Tafalla 

M., 2001, A&A, 372, 899 
Beltran M.T., Gueth F., Guilloteau S., Dutrey A., 2004, 

A&A, 416, 631 
Benedettini M., Yates J.A., Viti S., Codella C, 2006, MN- 

RAS, 370, 229 

Cesaroni R., Neri R., Olmi L., Testi L., Walmsley CM., 

Hofner P., 2005, A&A, 434, 1039 
Charnley S.B., Tielens A.G.G.M., Millar T.J., 1992, ApJ, 

399, L71 

Collings M.P., Anderson M.A., Chen R., Dever J.W., Viti 
S., Williams D., McCoustra M.R.S., 2004, MNRAS, 354, 
1133 

Danby G., Flower D.R., Valiron P., Schilke P., Walmsley 

CM., 1988, MNRAS, 235, 229 
Davis C.J., Eisloffel J., 1995, A&A, 330, 851 
Gueth F., Guilloteau S., Bachiller R., 1996, A&A, 307, 891 
Gueth F., Guilloteau S., Dutrey A., Bachiller R., 1997, 

A&A, 323, 943 
Gueth F., Guilloteau S., Bachiller R., 1998, A&A, 333, 287 
Hodapp K.-W., 1994, ApJSS, 94, 615 

Jimenez-Serra I., Martin-Pintado J., Rodriguez-Franco A., 

Marcellino N., 2004, ApJ, 603, L49 
Nisini B., Codella C, Giannini T., Santiago Garcia J., 

Richer J.S., Bachiller R., Tafalla M., 2007, A&A, 462, 163 
Perez-Gutierrez, M., 1999, Ph. D. "Quimica de choques en 

el flujo bipolar L1157", Universidad Complutense, Madrid 
Pineau des Forets G., Flower D.R., Chieze J. P., 1997, in 

Reipurt B., Bertoux C, eds, Proc. IAU Symposium n. 

182, Herbig-Haro flows and the birth of stars, p. 199 
Schilke P., Walmsley CM., Pineau des Forets C, Flower 

D.R., 1997, A&A, 321, 293 
Tafalla M., Bachiller R., 1995, ApJ, 443, L40 
Umemoto T., Iwata T., Fukui Y., Mikami H., Yamamoto 

S., Kameya O., Hirano N., 1992, ApJ, 392, L83 
Viti S., Codella C, Benedettini M., Bachiller R., 2004, MN- 
RAS, 350, 1029 
Zhang Q., Ho P.T.P., Wright M.C.H., Wilner D.J., 1995, 

ApJ, 451, L71 

Zhang Q., Ho P.T.P., Wright M.C.H., 2000, ApJ, 119, 1345 



